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Abstract 



We consider Kirchhoff equations with a small parameter e in front of the second-order 
time-derivative, and a dissipative term whose coefficient may tend to as i — > -|-oo 
(weak dissipation). 

In this note wc present some recent results concerning existence of global solutions, 
and their asymptotic behavior both as t ^ +00 and as e 0"^. Since the limit 
equation is of parabolic type, this is usually referred to as a hyperbolic-parabolic singular 
perturbation problem. 

We show in particular that the equation exhibits hyperbolic or parabolic behavior 
depending on the values of the parameters. 
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1 Introduction 



Let if be a separable real Hilbert space. For every x and y in H, \x\ denotes the norm 
of X, and {x, y) denotes the scalar product of x and y. Let A be a self-adjoint linear 
operator on H with dense domain D{A). We assume that A is nonnegative, namely 
{Ax,x) > for every x e D{A), so that for every a > the power A'^x is defined 
provided that x lies in a suitable domain D{A"). 

Let b : [0, +oo) — > (0, +oo) and m : [0, +oo) — > [0, +oo) be two given functions. For 
every £ > we consider the Cauchy problem 

+ b{t)u',{t) + m{\A^/\e{t)\^)Aue{t) = 0, (1.1) 

u{0) = uo, u'{0) = Ml. (1.2) 

This is the dissipative version of the celebrated equation introduced by G. Kirchhoff 
in [19] as a simplified model for transversal vibrations of elastic strings. We refer to the 
survey [14] for the non-dissipative case where e — 1 and b{t) = 0. Let us set 

li:= inim(a), 5:= inf b(t), v inf | : x g D(A), x ^ 

o->0 t>0 [ \xY 

Several features of (1.1) depend on the values of /x, 5, v. Let us recall some standard 
terminology. 

• 1^ on- degenerate vs degenerate equations These terms refer to the nonlinearity. 
Equation (1.1) is called nondegenerate or strictly hyperbolic if /i > 0, and degener- 
ate or weakly hyperbolic if > 0. The Cauchy problem (1.1), (1.2) is called mildly 
degenerate if > but 

m{\A'/\o\') ^ 0. (1.3) 

Whenever we consider degenerate equations, we always limit ourselves to the 
mildly degenerate case. The really degenerate case where m(\A^I'^UQ\^^ — seems 
to be still quite unexplored. 

• Constant vs weak dissipation We have constant dissipation when b{t) — 5 > for 

every t > 0, and weak dissipation when b{t) — > as t ^ +oo. Almost all known 
results for the constant dissipation case can be easily extended to non-constant 
dissipation coefficients provided that 6 > and b'{t) is bounded. For simplicity 
we often limit ourselves to the model case where b{t) = {1 + t)~P for some p > 0, 
the case p — corresponding to constant dissipation. 

In this note we don't consider equations with strong dissipation, which usually 
refers to dissipative terms of the form A°'u'^{t) with o; > 0, or better a > 1/2. 

• Coercive vs non-coercive operators The operator A is called coercive when u > 0, 
and it is called noncoercive when u > 0. This property of the operator has a great 
influence on the asymptotic behavior of solutions. 
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The singular perturbation problem in its generality consists in proving the conver- 
gence of solutions of (1.1), (1-2) to solutions of the first order problem 

b(t)u'^(t) + m{\A^/\,{t)f)Aue(t) = 0, u{0) = uo, (1.4) 

obtained setting formally e = in (1.1), and omitting the second initial condition 
in (1.2). Following the approach introduced by J. L. LlONS [20] in the linear case, one 
defines the corrector 9e{t) as the solution of the second order linear problem 

eO'^it) + b{t)e'^{t) =0 Vt>0, 

^,(0) = 0, ^^(0) = ^1 + -^m{\A'/^Uo\'')Auo w^. (1.5) 

It is easy to see that ^^^(0) = ^^(0) — u'{0), hence this corrector keeps into account 
the boundary layer due to the loss of one initial condition. Finally one defines r^it) and 
Ps{t) in such a way that 

Ue{t) = u{t) + 9,{t) + r,{t) = u{t) + pS) > 0. (1.6) 

With these notations, the singular perturbation problem consists in proving that 
r^{t) — > or p^{t) — > in some sense as £ — > O"*". The general problem can be spht into 
at least six subproblems. 

(1) Parabolic problem: global existence and decay estimates This is the first and 
usually easiest step in the theory. It consists in proving that (1.4) admits a unique 
global solution u{t), and then in estimating its decay rate as t — +oo. This 
decay rate is afterwards used as a benchmark for the decay rate of solutions of the 
hyperbolic problem. 

(2) Local existence for the hyperbolic problem and local-in-time error estimates Let 
T > be fixed. This subproblem consists in proving that, for every e > small 
enough, the solution Me(t) of the hyperbolic problem (1.1), (1.2) is defined (at least) 
on the interval [0,T], and in this interval Ue{t) converges to the solution u{t) of 
the limit problem. In this case the smallness of £, as well as the convergence rates, 
may depend on T. 

(3) Hyperbolic problem: global existence This subproblem consists in proving that 
problem (1.1), (1.2) admits a global- in-time solution provided that e > is small 
enough. From the point of view of existence, this is a strengthening of the previous 
step, and of course in general it requires stronger assumptions. 

Existence of global solutions without the smallness assumption on e is a widely 
open question, which seems to be as difficult as the non-dissipative case (see Sec- 
tion 4). 
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(4) Hyperbolic problem: decay estimates Once we know that the hyperbolic problem 
admits a global-in-time solution Us{t), a natural question concerns its behavior as 
t — > -|-oo {e is now small and fixed). What one expects in reasonable situations 
is that Us{t) decays as the solution u{t) of the corresponding parabolic equation. 
This is what has been actually proved in many cases. 

(5) Singular perturbation problem: global-in-time error estimates This subproblem 
is just the global-in-timc version of subproblem (2). The goal is therefore to give 
time-independent estimates on ps{t) or r^lt) as £ — > 0"*". 

(6) Singular perturbation problem: decay- error estimates This is the meeting point 
of subproblems (4) and (5), and it is the ultimate goal of the theory. It consists 
in estimating in the same time the behavior of Ui;{t) as t — -|-oo and as £ — >■ O"*". 
The general form of a decay-error estimate is something like 

<^(£)7(t) or |yl"r,(t)| <^(£)7(t). 

Of course one expects jit) to be the decay rate of solutions of the parabolic 
problem (or even better), and uj[e) to be the convergence rate which appears in 
the local-in-time error estimates. 

This program has generated a considerable literature in the last thirty years, for 
which we refer to the introductions of the following sections. In this note we sum up 
the state of the art and the main open questions. A rough overview is provided by 
Table 1, where we show, under different assumptions, which subproblems have received 
a reasonable or partial answer up to now. We focus in particular on the model dissipation 
coefficient of the form b{t) := (1 -|- t)~^ with p > 0, and on nonlinear terms which are 
either non-degenerate or of the form m{a) = a'^ for some 7 > (note that we allow also 
the non-Lipschitz case 7 e (0, 1)). 





p = 


p e (0. 1]. // > 


p e (0. 1]. > 


At > 


1-2-3-4-5-6 


l-2-3-4p-5-6p 


1-2-3-4-5-6 


m((j) = a'', 7 > 1 


l-2-3-4-5p 


1-2-3-4 


l-2-3p-4p 


m{a) = a^,-f e (0, 1) 


1-2-3-4 


1-2-3-4 


l-2-3p-4p 


m(cr) > (Lip. cont.) 


1-2-3-4 


1-2 


1-2 



Table 1: state of the art on subproblems (1) through (6). Numbers refer subproblems, 
a "p" means that in that case only partial (non-optimal) results have been obtained 

Looking at Table 1, one can guess that p = 1 plays a special role in the theory. This 
is true also in the linear case. Let us indeed consider equation 

au"{t) + --^—u'it) + cAuit) = 0, (1.7) 
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where a, b, c are positive parameters, and p > 0. This equation was investigated by 
T. Yamazaki [32] and J. Wirth [30]. They proved that (1.7) has both parabohc and 
hyperbohc features, and which nature prevails depends on p. When p < 1 the equation 
has parabohc behavior, in the sense that all its solutions decay to as i — > +oo as 
solutions of the parabolic equation with a = 0. When p > 1 the same equation has 
hyperbolic behavior, meaning that every solution is asymptotic to a suitable solution of 
the non-dissipative equation with 6 = (and in particular all non-zero solutions do not 
decay to zero). In the critical case p — 1 the nature of the problem depends on b/a, 
with the parabolic behavior prevailing as soon as the ratio is large enough. 

Our results for Kirchhoff equations are consistent with the linear theory. Indeed we 
have always hyperbolic behavior when p > 1, meaning that non-zero global solutions 
(provided that they exist) cannot decay to 0. Whenp < 1 we were able to prove that the 
behavior is of parabolic type in many cases. In all such situations the critical exponent 
p — 1 falls in the parabolic regime, but this is simply due to the fact that in our equation 
we have that 6 = 1 and a = e is small enough, hence the ratio b/a is always big enough. 

For shortness's sake we don't include proofs in this note. Nevertheless, we conclude 
this introduction by mentioning the useful energies and the technical reasons why the 
problem becomes harder and harder when the equation is degenerate, the dissipation is 
weak, and the operator is non-coercive. In the parabolic case all estimates follow from 
the monotonicity of the classical energies 

E,it):^\A'/Mt)\', Pit):- 



In the hyperbolic case, all known techniques for proving global existence for (1.1) 
require an a priori estimate such as 

If // > and b{t) is a positive constant, an a priori bound on Aue{t) and u'^{t), 
together with the smallness of £, is enough to establish (1.8). When b{t) ^ as 
t +00, the boundedness is no more enough, and we need some a priori informations 
on the decay of AUfr{t) and u'^{t). This means that global existence and decay estimates 
become intimately tied, and they must be treated together. The main energies involved 
in these estimates are 

E,,{t) := e^l^^^ + \A('^'y\M^ G^t) := (^ g) 

where ^(t) := m{\A^/'^Ue{t)\'^). They are both extensions of the first energy of the 
parabolic case. The use of E^_k{t) is quite classical, and dates back to [3, 31], while 
Gs{t) was introduced by the authors in [8]. 
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The degenerate case is more complex. Let us assume for example that 'm{a) = if 
and only if cr = 0. Then the decay of the solution implies that the denominator in the 
left-hand side of (1.8) tends to 0. When m{a) — with 7 e (0, 1), then also the term 
with vnl diverges to +00 as the solution approaches 0. This comphcates proofs both in 
the case of constant, and in the case of weak dissipation. 

The basic idea to deal with degenerate nonlinear terms of the form m(a) = cr'"' is to 
exploit that am\a)/m{a) has a finite limit as a — > O"*". This reduces (1.8) to 



\AuM-WM 



<h{t). (1.10) 







to 






to| 



This inequality has been approached using (1.9) and the further energies 

Ce |AV2^e|^ |AV2m^|2' 

These energies have been introduced by the first author in [7] as a hyperbolic version 
of the second energy of the parabolic case. 

If m{a) is a general nonnegative (even Lipschitz continuous) function, then it may 
happen that am' {a) / m{a) is unbounded in a neighborhood of cr = 0. In this case (1-10) 
does not imply (1.8), and each step seems to require new ideas. For this reason we are 
quite skeptic about a future relevant progress in the last line of Table 1. 

Concerning the coerciveness of the operator, it is well known that small eigenvalues 
deteriorate the decay estimates on solutions, even in the parabolic case, and we have 
seen that decay estimates are fundamental also for existence issues. 

This note is organized as follows. Section 2 is devoted to subproblem (1), namely 
existence and decay estimates for the parabolic problem. Section 3 is devoted to sub- 
problem (2), namely local existence results for the hyperbolic equation and local- in-time 
error estimates for the singular perturbation problem. In Section 4 we show that (1.1) 
has hyperbolic behavior whenever p > 1. Section 5 is devoted to global existence and 
decay-error estimates in the nondegenerate case. Section 6 is devoted to the degenerate 
case. Finally, Section 7 is a collection of open problems. 



2 The parabolic problem 

The theory of parabolic equations of Kirchhoff type is quite well established. This 
equation appeared for the first time in the pioneering paper [2] by S. Bernstein. He 
considered the concrete equation in the interval (0, 1), with a nondegenerate nonlinearity 
and constant dissipation, and he proved that for every initial condition in the Sobolev 
space i/^((0, 1)) the equation admits a unique solution, which is actually analytic in the 
space variable for every t > (the classical regularizing effect of parabolic equations). 
This result was afterwards extended by many authors (see [1, 21]). 
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The more general version is probably stated in [15]. The basic fact observed in [15] 
is that any solution u{t) of (1.4) can be written in the form 

u{t) — v{a{t)), 

where v{t) is the solution of the linear Cauchy problem with constant coefficients 

v'(t)+Av{t)^0, v(0)^uo, (2.1) 

and a : [0,+oo) — > [0, +cx)) is the solution of the ordinary differential equation 

b{t)a'{t) = m{\A^/\{a{t))\^), a{0) = 0. 

In other words, the solution of (1.4) is always a time reparametrization of the solution 
of the heat-like equation (2.1). At this point it is quite easy to prove the following 
existence result. 

Theorem 2.1 (Global existence for the parabolic problem) Let H be a Hilbert 
space, and let A be a nonnegative self-adjoint (unbounded) operator on H with dense 
domain. Let m : [0, +00) — > [0, +00) be a locally Lipschitz continuous function, and let 
uoEDiA). 

Then problem (1-4) has a unique global solution 

ueC^ ([0, +00); H) n C° ([0, +00); D{A)) . 

If in addition m{\A^/'^uo\'^)Auo 7^ 0, hence u'{0) ^ 0, then the solution is non- 
stationary, and u & ((0, +00); Di^A")) for every a >0. 

Decay estimates for u{t) can be deduced from decay estimates for (2.1) and the 
asymptotic behavior of the parametrization a(t). Concerning (2.1), it is well known 
that the asymptotic behavior of solutions depends on the coerciveness of the operator 
A. If A is coercive with some constant v > 0, then solutions decay exponentially to 0, 
with a rate depending on u. In this case we have indeed that 

\A'/'uo\'e^p (^-2^^^t) < \A'/Mt)\' < \A'/'uo\'eM-'^^t). 

If A is non-coercive {v > 0), then decay rates are slower. We have indeed that 



Note in particular that the estimates from below and from above for |74^/2t'(t)|2 
involve different rates. This range of rates cannot be improved because, when the 
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i/ > 




o 

A 




ci(l + tfPe-^(^+'r+' < < C2(l + i)2Pe-"^(i+*)''+^ 


Al 


i/ > 






iy>0 


(l+^)(p+l)/7 ^ 1^"^ ^ (l+i)(p+l)/7 
(l + t)(p+l)/7 ^ I^^^^Wl ^ (l + ^^(p+l)/7 


<- 

K 




(1 + t)2+(p+l)/7 - \^^^>\ - (1 + t)2+(p+l)/7 




(1 + t)(p+l)/7 ^ 1 "^^''1 ^ (1 + 






\Au{t)\^ < , 

1 Wl - (1 _^ i)(P+l)/7 
|m {t)\ ^ ^^[2^2+(i_p)^+p+i)/(^2+^) 



Table 2: Decay estimates for the parabolic problem 



operator has a sequence of eigenvalues converging to 0, any intermediate rate is realized 
by a suitable solution. 

Once we know the decay of v{t), we can easily deduce the asymptotic behavior of 
a{t), hence also the asymptotic behavior of u{t). In Table 2 we sum up the decay esti- 
mates which can be obtained in this way, limiting ourselves for simplicity to dissipation 
coefficients of the form b{t) = (1 + t)~^ with p > 0, and to nonhnear terms which are 
either non-degenerate or of the form m{a) = a'^ with 7 > 0. 

We stress that in all these cases solutions decay to zero, and the decay rate be- 
comes stronger and stronger as p grows. This contrasts with the hyperbolic case, where 
solutions cannot decay when p > 1 (sec Section 4). 
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3 Local-in-time error estimates 



All the local existence results for the non-dissipative equation (see [14, Theorem 2.1]) can 
be easily extended to the dissipative case. This provides a continuum of local existence 
results, with the regularity requirements on the initial data depending on the continuity 
modulus of m. In this note we limit ourselves to Lipschitz continuous nonlinear terms, 
or to the non- Lipschitz case m{a) = with 7 e (0,1), where the nondegeneracy 
assumption (1.3) makes the problem just mildly non- Lipschitz. In all these cases the 
equation is locally well posed for initial data in Sobolev spaces. 

In this section we focus on a property which is slightly stronger than local existence, 
and which could be called almost global existence. The first result is indeed that the life 
span of Ue(t) tends to -|-cx) as £ ^ O"*". 

Theorem 3.1 (Hyperbolic problem: almost global existence) Let H be a Hil- 
bert space, and let A be a nonnegative self-adjoint (unbounded) operator on H with 
dense domain. Let m : [0, +00) — > [0, +00) and b : [0, +00) —>■ (0, +00) be two locally 
Lipschitz continuous functions. Let us assume that {uo,Ui) e D{A) x D{A^/'^) satisfy 
the non- degeneracy assumption (1.3), and let T > 0. 

Then there exists £0 > such that for every e e (0,£o) problem (1.1), (1-2) has a 
unique solution 



Then we study the convergence of Ueit) to the solution u{t) of the limit problem. 

Theorem 3.2 (Singular perturbation: local-in-time error estimates) Let H , 

A, m{a), b{t), uq, ui, T, Eq be as in Theorem 3.1. Let u{t) be the solution of the 
corresponding parabolic problem (I.4), o,nd let rs{t) and Ps{t) be defined by (1.6). 
Then we have the following conclusions. 

(1) Without further assumptions on initial data, hence {uo,ui) G D{A) x D^A^^"^), we 



(2) If in addition we assume that {uo,Ui) e D{A'^^'^) x D{A^^'^), then there exists a 
constant C such that for every e e (0, £0) we have that 



u, e CmO,T];H)nC^ {[0,T];D{A^/^))nC°{[0,T];D{A)). 



have that 



p,{t)\^ + \A'/^Pe{t)\^ + \ApM^ + \Kit)\^ ^ uniformly m [0,T], 




pM' + \A'/'Pe{t)\' + e\r'M' < Ce' Vt e [0,T], 




8 



(3) If in addition we assume that {uo,ui) G D^A^) xD{A), then there exists a constant 
C such that for every e G (0, eo) we have that 



\ApM' + \r'M' + e\A'/'r'M' < Ce' G [0,r], 

r \A'/VM'' dt < Ce\ 
Jo 

We point out that the remainder r^{t) is well suited for estimates involving deriva- 
tives, because it doesn't feel the effects of the boundary layer due to the loss of one initial 
condition. On the contrary, the remainder Ps{t) is better suited for estimates without 
derivatives. This is because, for example, Ap^{0) is defined whenever uq G D[A), while 
Are{0) requires uq G D{A'^) (see definition (1.5) of wq). 

Both the existence and the convergence result are local-in-time, namely constants, 
error estimates, and the smallness of e do depend on the interval [0, T] chosen at the 
beginning. On the other hand, the assumptions required on b{t) and m{(T) are quite 
weak. The dichotomy between hyperbolic and parabolic behavior mentioned in the 
introduction appears only as t — > -|-oo, hence it plays no role on a fixed time interval. 
In particular we don't need to assume that p < 1 in the case where b{t) = (1 + 

Prom Theorem 3.2 it is clear that convergence rates for the singular perturbation 
problem depend on the regularity of initial data. This situation is consistent with the 
linear case. Indeed in [9] we considered the linear equation 

eu"{t) + u'^{t) + Aue{t) = 

and the corresponding limit parabolic problem, and we proved similar results. We 
also proved that an error estimate such as \A^/'^ pe{t)\'^ < Ce^ is possible only when 
(mo,mi) G D{A^I^) X D{A^'^). 

We just remark that uniform convergence, without any rate, requires initial data 
in spaces such as D{A) x D^A^^"^), hence with "gap 1/2" between the regularity of 
u'^ and Ue, a typical feature of hyperbolic problems. On the contrary, if we want some 
convergence rate, we have to work in spaces such as D{A^/'^) xD{A^/'^) oi D{A^)xD [A], 
hence with "gap 1", a typical feature of parabolic problems. In our opinion, this gives 
further evidence that the parabolic nature dominates in the limit. 

A formal proof of Theorem 3.1 and Theorem 3.2, as they are stated, has never 
been put into writing. Error estimates have been considered in at least three papers, 
always with constant dissipation. B. F. Esham, and R. J. Weinacht [5] proved error 
estimates in the nondegenerate case with initial data in D{A^/'^) x D{A). The second 
author [16] considered the degenerate case, proving uniform convergence for data in 
D{A) X D(Ai/2), and error estimates for more regular data. Finally, the authors [11] 
proved error estimates in the degenerate case with the optimal requirement that initial 
data are in D(A^/^) x D(A^/^) (see [11, Proposition A.lj). It should be quite standard to 
extend those proofs to the case of weak dissipation, just because on a fixed time interval 
the function h{t) is always strictly positive. 
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4 The hyperbolic regime 



In this section we show that when the dissipation is too weak, namely b{t) too 
fast, then equation (1.1) behaves in a hyperbohc way, in the sense that its non-zero 
global solutions (provided that they exist) do not decay to as t ^ +oo. Of course 
this doesn't prevent such solutions from existing (which remains an open problem), but 
it shows that the problem cannot be approached using the standard methods based on 
estimates such as (1.8) or (1.10). Since solutions of the parabohc problem always decay 
to 0, this shows also that no decay-error estimate can be true in this case. Note that 
condition (4.1) is equivalent to p > 1 when b{t) = (1 + i)~^. 

Theorem 4.1 (Hyperbolic regime) Let H be a Hilbert space, and let A be a non- 
negative self-adjoint (unbounded) operator on H with dense domain. 

Let m : [0, +oo) [0, +oo) he a continuous function. Let b : [0, -|-oo) (0, -|-oo) be 
a continuous function such that 

r+oo 

I b(s) ds < +00. (4.1) 
Jo 

Let {uo,Ui) e D{A) x D{A^/^) be such that 

\ui\^+ / m{a)da > 0. (4.2) 



Let us assume that for some e > problem (1.1), (1-2) has a global solution 

u, e C'([0, +oo); H) n C^([0, +oo); D{A^'^)) n C°([0, +oo); D{A)). (4.3) 
Then we have that 

liminf (|<(t)|' + |A^/V(t)r) > 0. (4.4) 

t— >+oo 

The proof of this result is very simple, and relies on the usual Hamiltonian 

»|AV2„,(t)|2 



Heit) -.^ e\u'^{t)\'^ + m{a)da. 
Jo 



Assumption (4.2) is equivalent to say that 7i£{0) > 0. Moreover we have that 

n',{t) = -2b{t)\u'M'' > ~b{t)ne{t) yt > o, 

hence ^ 

ne{t) > HeiO) exp ( -- [ b{s) ds] \/t > 0. 



For a fixed £ > 0, the right-hand side is greater than a positive constant independent 
on t because of (4.1) and the fact that He{0) > 0. This imphes (4.4). 
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5 The nondegenerate case 



In this section we focus on the hyperbohc equation (1.1) under the non-degeneracy 
assumption fj, > 0. 

The case with constant dissipation was considered independently by E. H. de Brito [3] 
and by Y. Yamada [31]. They proved existence of a global solution provided that e is 
small enough. Decay estimates for these solutions were proved by Y. Yamada [31] in 
the non-coercive case, and by E. H. de Brito [4] and by M. Hosoya and Y. Yamada [18] 
in the coercive case. All these estimates were afterwards reobtained as a particular case 
of the theory developed in [10]. 

More recently, H. Hashimoto and T. Yamazaki [17] proved that for initial data 
(mo,mi) e X D{A) one has that 

\p,{t)\^ + (1 + t)\A'/^pM^ + e{l + t)yM^ < Ce" > 0, 

where of course C doesn't depend on t and e. When {uq,Ui) G D{A^) x D{A), the 
coefficient e in the left-hand side may be dropped, thus providing a better convergence 
rate on r'^{f^. This is a first example of decay-error estimate. 

The weakly dissipative case was considered only in last years. Apart from a result 
obtained in a special situation by M. Nakao and J. Bae [24], the problem in its full 
generahty was solved by T. Yamazaki [33] in the subcritical case p < 1 with some 
technical requirements on initial data, and then by the authors [12] (see also [34]) in the 
general case p <1 with minimal requirements on initial data. 

The results are the following. 

Theorem 5.1 (Hyperbolic problem: global existence) Let H he a Hilbert space, 
and let A be a nonnegative self-adjoint (unbounded) operator on H with dense domain. 
Let fi > 0, and let m : [0, +oo) [/i, +oo) be a locally Lipschitz continuous function. 
Let hit) := {l+t)-P with p e [0,1], and let (mo,mi) e D{A) x D{A^/^). 

Then there exists £o > such that for every e e (0, £o) problem (1.1), (1-2) has a 
unique global solution satisfying (4-3). 

Theorem 5.2 (Hyperbolic problem: decay estimates) Under the same assump- 
tions of Theorem 5.1 there exists a constant C such that for every e e (0,£o) have 
that 

\uM^ + (1 + tY+'\A'/\,{t)f + (1 + tfiu'Mf <c yt>o, 

e\A'/\'M' + \Aus{t)\' < ^ yt > 0, 

r+oo 

/ {i + ty{\u'M' + \A'^Wt)n dt<c, 

Jo 

r+oo 

/ (1 + t)2f+i (|AVX(t)|2 + \AuM^) dt < c. 
Jo 
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Theorem 5.3 (Singular perturbation: decay-error estimates) Let H , A, fi, 

m{a), b{t), p, uq, ui, £q he as in Theorem 5.1. Let u(t) be the solution of the corre- 
sponding parabolic problem (1-4), CLnd let re{t) and pe{t) be defined by (1.6). 
Then we have the following conclusions. 

(1) Without further assumptions on initial data, namely {uq.,ui) e D[A) x D^A^/"^), 
we have that 

\pM? + (1 + trV^VWI' + (1 + tf'^+^~^\ApS)? + (1 + tf\r'M? ^ 
uniformly in [0,+oo), and 



I 



f 

Jo 

+00 



{i + ty{\r'M' + \A'^'pemdt^o, 
(1 + tfp+' (|A'/V;(OI' + dt ^ 0. 



'0 

(2) If in addition we assume that {uq,ui) G D{A^^'^) x D{A^/^), then there exists a 
constant C such that for every e G (0, £o) we have that 

+ (l + i)f+V'''VeWr + £(l + ir+>eWr < Ce" \ft > 0, 

r+oo 

/ {l + tr i\r'M' + \A'^'Pe{mdt <Cs\ 
Jo 

(3) If in addition we assume that (mq, ui) e D{Al^) xD{A), then there exists a constant 
C such that for every e e (0, £o) we have that 

(1 + t)2(f+i)|Ap.Wr + (1 + tfWeif)? < Ce'' \/t > 0, 



/" 

Jo 



(1 + tfP+' {\A^/^r',{t)f + \Ap,{t)f) dt < Ce\ 

lo 

We point out that in the previous three theorems the operator A is never assumed 
to be coercive. 

The global-in-time convergence rates (with respect to e) appearing in Theorem 5.3 
are optimal because they coincide with the local-in-time convergence rates of Theo- 
rem 3.2, which in turn are the same of the linear case. 

The decay rates (with respect to time) appearing in Theorem 5.2 and Theorem 5.3 
are optimal for non-coercive operators. In this case indeed they coincide with the decay 
rates of solutions of the corresponding parabolic equation, as shown in Table 2. 

In the coercive case these decay rates are not optimal. In the case p = indeed we 
know that solutions exponentially decay to zero as solutions of the parabolic problem 
(see [4, 18, 10]). We strongly suspect that the same is true also for every p e [0,1], 
namely that solutions decay as shown in the first three rows of Table 2. Of course 
also the decay rates in Theorem 5.3 should be changed accordingly. We give no precise 
statement or reference because this part of the theory has never been put into writing. 
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6 The degenerate case 



Several papers have been devoted to global existence and decay estimates for equation 
(1.1) in the degenerate case // > 0. Let us begin with constant dissipation. In this case 
global existence results (provided that the problem is mildly degenerate and e is small 
enough) were proved by K. Nishihara and Y. Yamada [25] in the case where m{a) = a'^ 
(with 7 > 1), by the authors [8] in the case where m{a) > is any Lipschitz continuous 
function, and by the first author [6, 7] in the non-Lipschitz case where m((T) = a'^ with 
7 6(0,1). 

Decay estimates have long been studied for equations with constant dissipation. 
In the case m((T) = a'^ with 7 > 1, the first decay estimates were obtained by K. 
Nishihara and Y. Yamada [25] in the coercive case, and by K. Ono [29] in the non- 
coercive case. The case m{a) — with 7 e (0,1) was considered in [6]. In the 
special case m((T) = g, T. Mizumachi [22, 23] and K. Ono [26, 27] proved better decay 
estimates, namely estimates with decay rates which are faster than those obtained by 
putting 7 = 1 in the previous ones. This in particular showed that the previous results 
were not optimal. 

A complete answer was given by the authors in [10], where the case of a general 
nonlinearity m{a^ > is considered. The decay rates obtained in [10] coincide with the 
decay rates of solutions of the parabolic problem. 

Let us consider now the equation with weak dissipation, focussing on the model case 

£<(t) + ^^uXt) + \A'l^uM^ Au,{t) = 0, (6.1) 

of course with the mild non-degeneracy assumption (1.3). The only previous result we 
are aware of was obtained by K. Ono [28]. In the special case 7 = 1 he proved that a 
global solution exists provided that e is small and p E [0, 1/3]. The reason of the slow 
progress in this field is hardly surprising. In the weakly dissipative case existence and 
decay estimates must be proved in the same time. The better are the decay estimates, 
the stronger is the existence result. 

Ten years ago decay estimates for degenerate equations were far from being optimal, 
but for the special case 7 = 1. In [10] a new method for obtaining optimal decay 
estimates was introduced. This allowed a substantial progress on equation (6.1). 

Let us begin with our existence and decay results proved in [13]. The first one 
concerns the coercive case. 

Theorem 6.1 (Coercive case: global existence and decay estimates) Let H be 

a Hubert space, and let A be a nonnegative self-adjoint (unbounded) operator with dense 
domain. Let us assume that A is coercive (u > 0). Let 7 > 0, and let p G [0, 1]. Let us 
assume that {uo,ui) e D{A) x D{A^/^) satisfy (1.3). 

Then there exists > such that for every e e (0,£o) problem (6.1), (1.2) has a 
unique global solution satisfying (4-3). 
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Moreover there exist positive constants Ci and C2 such that 



c c 

|2 / ^2 



\uM\ < ^-^ ^ ^^2+{p+l)/7 



> 0. 



We point out that Theorem 6.1 is optimal both in the sense that all p G [0, 1] are 
considered, and in the sense that solutions decay as in the parabohc case (see Table 2). 
In the non-coercive case we have the following result. 

Theorem 6.2 (Non-coercive case: global existence and decay estimates) Let 

H he a Hilbert space, and let A be a nonnegative self-adjoint (unbounded) operator with 
dense domain. Let 7 > 1; and let 

7^ + 1 
72 + 27 - 1 ■ 

Let us assume that {uo,ui) e D{A) x D{A^/^) satisfy (1.3). 

Then there exists £0 > such that for every e e (0,£o) problem (6.1), (1.2) has a 
unique global solution satisfying (4-3). 

Moreover there exist constants Ci and C2 such that 



0<P< 2,0.. (6-2) 



(1 + t)(p+i)/7 - I ^1 ~ (1 + t)(p+'^y('y+^) 
\AuJt)\'^ < ^f-TT^ > 0, 

I e\ /I - (^IJ^ ^^(p+l)/7 - ' 

l^eWI ^ (1 + i)[272+(l-p)7+P+l]/(72+7) - 

Theorem 6.2 doesn't represent a final answer in the non-coercive case. Let indeed 
denote the right-hand side of (6.2). It is easy to see that p-y < 1 for every 7 > 1, with 
equality only when 7 = 1, and asymptotically as 7 ^ +00. Since we have hyperbolic 
behavior when p > 1 (see Section 4), and parabolic behavior for p G [0,j9^], this means 
that there is a non-man's land between p^ and 1 where things are not clear yet. 

The only case where this region is empty is when 7 = 1. In this case all exponents 
p e [0, 1] fall in the parabolic regime, and this improves the result obtained in [28] 
(p G [0, 1/3]) also in the case m{(j) = a. 

We stated Theorem 6.2 assuming 7 > 1. In the case 7 G (0, 1) we have a weaker 
result, namely global existence for p G [0,7/(7 + 2)] (see [13, Remark 2.6]). Figure 1 
represents hyperbolic and parabolic regimes, and the no-man's land in between. 

The singular perturbation problem is still quite open in the degenerate case. We 
have indeed only the following partial result for the constant dissipation case (see [11]). 
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Figure 1: parabolic and hyperbolic regimes in the degenerate non-coercive case 



Theorem 6.3 (Constant dissipation: global-in-time error estimates) Let H he 

a Hilbert space, and let A be a nonnegative self-adjoint (unbounded) operator with dense 
domain. Let Ugit) be the solution of equation (6.1) with 7 > 0, p = 0, and initial data 
{uq.ui) E D(A) X D^A^^"^) satisfying (1.3). Let u(t) be the solution of the corresponding 
parabolic problem, and let r^it) and peit) be defined by (1.6). 
Then we have the following conclusions. 

(1) Without further assumptions on initial data, namely {uo,ui) E D{A) x D{A^^'^), 
we have that 

\pM' + \A'^'Peit)\' + \ApM' + Kit)]' 
uniformly in [0,+oo), and 

r+co 

Jo 

(2) If in addition we assume that 7 > 1 and (mo,mi) E x D{A^/'^), then there 
exists a constant C such that for every e E (0,£o) we have that 

\pe{t)\^ + e\A^'^Pe{t)\''<Ce'' Vt>0, 

/ \r'Xt)\^dt<Ce. 
Jo 

(3) If in addition we assume that 7 > 1 and {uo,ui) E D^A^^"^) x D{A), then there 
exists a constant C such that for every e E (0,eo) we have that 

Theorem 6.3 is far from being optimal. First of all most of the convergence rates in 
the second and third statement are weaker than the corresponding rates in Theorem 3.2. 
Moreover, all statements present just error estimates, and not decay-error estimates as 
in Theorem 5.3. It is possible to add some decays with some extra work, but in any 
case they are so far from those appearing in Theorem 6.1 and Theorem 6.2 that we 
decided not to include them. Last but not least. Theorem 6.3 is limited to equations 
with constant dissipation. 
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7 Open problems 



The main open problem in the theory of Kirchhoff equations is existence of global 
solutions. We have seen that in the dissipativc case an affirmative answer can be given 
provided that e is small enough. So the first question is whether this condition is 
necessary or not. 

Open problem 1 Let us consider equation (1.1) with m : [0, +cxd) [l,+oo) of class 
C°° , and constant dissipation b(t) = 1. Let us assume that {uq,Ui) G D{A°°) x D{A°°), 
where D{A^) is the intersection of all spaces D{A") with a > 0. 

Does the Cauchy problem (1-1), (1-2) admit a global solution for every £ > 0? 

We stated the question with generous assumptions both on the nonlinearity (smooth- 
ness and strict hyperbolicity), and on initial data (regularity). In any case there are no 
counterexamples, even with less regular terms and data, or with h{t) = 0. 

Even assuming the smallness of e, one may ask if a global solution exists under 
assumptions weaker than those required in the previous sections. This leads to the 
following question. 

Open problem 2 Let us consider the Cauchy problem (1-1), (1-2) in each of the fol- 
lowing situations. 

• In the hyperbolic regime where h{t) — (1 + with p > 1. 

• In the case where assumption b{t) = (l + t)"^ with p < 1 is replaced by the weaker 
condition (4-1) ■ 

• In the really degenerate case m{\A^/'^uo\'^) = 0. 

Is it possible to prove global existence provided that e is small enough? 

A third question related to global existence issues concerns the regularity of initial 
data. All existence results stated in the previous sections assume that {uq, Ui) G D{A) x 
D{A^^'^). On the other hand, the classical local existence results for the non-dissipative 
equation require the weaker assumption (mo,mi) G D{A^/^) x D{A^/^). Therefore a 
natural question is whether the global existence results for dissipative equations can be 
extended to this weaker class of data. 

In the constant dissipation case, it is not difficult to give an affirmative answer when 
/X > or when m{a) = with 7 > 2. On the contrary, the proof given in [8] for 
a general locally Lipschitz continuous non-linearity m{a) > seems to require in an 
essential way that (ito,iti) G D{A) x D(A^/'^). So the problem is the following. 

Open problem 3 Let us consider equation (1.1) with constant dissipation h{t) = 1, 
and with any locally Lipschitz continuous nonlinearity m{a) > 0. Let us assume that 
{uo,ui) G D(A^/^) X D{A^/^) satisfy the non- degeneracy condition (1.3). 

Does problem (1.1), (1-2) admit a global solution for every small enough e? 
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The last open question concerning existence is how to fill the no-man's zone left by 
Theorem 6.2 and described in Figure 1. 

Open problem 4 Let H , A, Uq, Ui he as in Theorem 6.2. Let us assume that either 
7 e (0,1) andp G (7/(7 + 2), 1], or that j > 1 andp e (^7,1], wherep^ is the right-hand 
side of (6.2). 

Does problem (6.1), (1-2) admit a global solution whenever e is small enough? 

All previous examples suggest that the answer should be affirmative, but a proof 
seems to require some new ideas. 

The singular perturbation problem is arguably the new frontier in this research field. 
This problem has been quite well understood only in the nondegenerate case, in which 
case, however, the decay rates are optimal only for non-coercive operators. A first open 
question is therefore the following. 

Open problem 5 Let the assumptions of Theorem 5.3 he satisfied. Let us assume also 
that the operator A is coercive (u > 0). 

Prove the same conclusions of Theorem 5.3 with all polynomial decay rates such as 
(1 + t)'^ replaced hy exponential decay rates of the form exp(Q;(l -|- ty~^^), where a is a 
suitahle constant. 

The singular perturbation problem is quite open in the degenerate case. One should 
try to extend Theorem 6.3 in order to allow weak dissipations, and involve better decay 
and convergence rates. An example of open question is the following. 

Open problem 6 Let H, A, 7, p, uq, ui, sq he as in Theorem 6.1. Let u{t) he the 
solution of the corresponding paraholic prohlem, and let re{t) and pe{t) he defined hy 
(1.6). 

Under the appropriate conditions on initial data, prove that there exists a constant 
C such that for every e e (0, Eq) we have that 

(1 + t)^+^^'^\A^l^pe{t)\'' < Ce^ yt > 0, 

(1 + ^)2+(p+l)/7|^/(^)|2 < ^^2 y^ > Q 

In this estimates we require on rs{t) and pe{t) the same decay rates {as t ^ +00) 
of u^lt) and u{t) separately, and we require the same convergence rates (as £ ^ 0+) of 
the local-in-time error estimates. We actually suspect that in the degenerate case the 
remainders rs{t) and p^{t) decay faster than Us{t) and u{t). 
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